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We apply spherical needlets to the Wilkinson Microwave Anisotropy Probe 5-year cosmic mi- 
crowave background (CMB) data_set, to search for imprints of nonisotropic features in the CMB 
sky. We use the needlets' localization properties to resolve peculiar features in the CMB sky and 
to study how these features contribute to the anisotropy power spectrum of the CMB. In addition 
to the now well-known "cold spot" of the CMB map in the southern hemisphere, we also find two 
hot spots at greater than 99% confidence level, again in the southern hemisphere and closer to the 
Galactic plane. While the cold spot contributes to the anisotropy power spectrum in the multipoles 
between ^ = 6 to ^ = 33, the hot spots are found to be dominating the anisotropy power in the 
range between i = 6 and ^ = 18. Masking both the cold and the two hot spots results in a reduction 
by about 15% in the amplitude of the angular power spectrum of CMB around i = 10. The result- 
ing changes to the cosmological parameters when the power spectrum is estimated masking these 
features (in addition to the WMAP team's KQ85 mask) are within the la errors published with 
the WMAP mask only. We also study the asymmetry between the angular power spectra evaluated 
on the northern and southern hemispheres. When the features detected by needlets are masked, we 
find that the difference in the power, measured in terms of the anisotropy variance between 1 = 4 
and ^ = 18, is reduced by a factor 2. We make available a mask related to needlet features for more 
detailed studies on asymmetries in the CMB anisotropy sky. 

PACS numbers: 98.80.-k, 98.80.Es 



I. INTRODUCTION 

Beyond the angular power spectrum of the cosmic 
microwave background (CMB) anisotropies, the high- 
sensitivity all-sky CMB maps produced by the Wilkinson 
Microwave Anisotropy Probe (WMAP)^ [1] have enabled 
detailed statistics studies to extract higher-order infor- 
mation. These studies include characterizations of asym- 
metries in the CMB sky [2] and the search for anomalous 
features in the anisotropy field. Previous analyses have 
shown evidences for alignments [3], asymmetry in the 
CMB statistics between northern and southern Galactic 
hemispheres [4], and features such as the "cold spot," a 
significant negative feature in the CMB map first identi- 
fied with wavelets [5] . 

Wavelets can be constructed both in real space on the 
sphere, for instance by means of the so-called stereo- 
graphic projection [6-9], and in harmonic space, for in- 
stance following the prescription described in [10, 11], see 
also [12]. There have been already several applications of 
wavelets to CMB data analysis, including tests for non- 
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Gaussianity and asymmetries [13-17], polarization anal- 
ysis [18], foreground subtraction [19], foreground compo- 
nent separation [20, 21], and point source detection in 
CMB anisotropy maps [22]. The reason for such a wide 
success can be motivated as follows: the comparison be- 
tween models and CMB data are primarily undertaken 
in the Fourier domain, where each multipole can be mea- 
sured separately with all-sky experiments. However, in 
the presence of sky cuts data analysis must take into 
account the coupling among multipoles induced by the 
mask; statistical studies consequently become much more 
challenging. Wavelets address these issues by providing 
the possibility to combine sharp localization properties 
in pixel space with a peaked window function in the har- 
monic space. 

In the latest two years, needlets have been proposed 
as a new wavelet system which enjoys several advantages 
over existing constructions. The introduction of needlets 
into the mathematical literature is due to [23], where 
their localization properties in pixel and harmonic space 
are established. The analysis of their statistical proper- 
ties for spherical random fields is first provided by [24], 
where uncorrelation is discussed and it is also shown how 
needlets can be used to implement estimators for the an- 
gular power spectrum or tests of non- Gaussianity. The 
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first application to WMAP data is due to [25], where 
needlets are used to estimate (cross-) angular power spec- 
tra in order to search for dark energy imprints on the cor- 
relation between large-scale structures and CMB. A full 
description of needlets and their potentiality for CMB 
data analysis is then provided by [26]; Ref. [27] inves- 
tigates the effect of different window functions in their 
constructions, whereas Ref. [28] provides further math- 
ematical results on their behavior for partially observed 
sky-maps. The previous results have been extended to 
angular power spectrum estimation in the presence of 
noise [29, 30], estimation of the bispectrum [31], fore- 
ground component separation [32] , and analysis of direc- 
tional data [33]; see also Refs. [34-36] for further devel- 
opments. 

Among the advantages of needlets compared to previ- 
ously adopted wavelets, we recall that they are compactly 
supported in harmonic space, so that they allow to focus 
the analysis on a specific set of multipoles, which can 
be completely controlled [23] . Needlet coefficients can be 
shown to be (approximately) uncorrelated both across 
different frequencies and, at high multipoles, at different 
locations in pixel space, which makes statistical analysis 
much more efficient; it is possible to provide an exact 
analytic expression for their variance and correlation, in 
terms of the underlying angular power spectrum [24, 25]. 
Needlets do not rely on any tangent plane approxima- 
tion, but they are directly embedded into the manifold 
structure of the sphere; they allow for direct reconstruc- 
tion formulas (a consequence of the fact that they make 
up a tight frame system) and they are computationally 
very convenient [26]. 

In this paper, we make use of needlets to further 
study features in the WMAP CMB maps. We focus on 
the large angular scales or, equivalently, on multipoles 
smaller than 200. We recover the cold spot that was pre- 
viously detected in WMAP data with wavelets, and we 
also detect other features, including two hot spots, which 
have so far received less attention. By masking these fea- 
tures, we study how the angular power spectrum of CMB 
anisotropics is modified. Given that these features are 
located in the southern hemisphere, we also discuss the 
extent to which these features could be responsible for 
the north-south asymmetry in WMAP data [2, 37, 38]. 
As is well known, this asymmetry has also drawn much 
interest in the theoretical community, since it could entail 
strong implications on the physical nature of primordial 
perturbations, including inflation [39]. By masking the 
low-^ features, we find that the difference in the CMB 
anisotropy variance between the two hemispheres is re- 
duced by a factor of 2, reducing the significance of pre- 
vious detections. 

We also explore the evidence that statistically signifi- 
cant bumps and dips in the CMB anisotropy power spec- 
trum at multipoles of 20 and 40 could be related to fea- 
tures in the CMB sky. We did not locate any partic- 
ular feature in the sky that dominates the power spec- 
trum at these multipoles; masking the significant features 



detected by needlets tuned to these multipoles did not 
change the power spectrum more than ~5%. 

The paper is organized as follows: In the following 
section we describe briefly the needlet formalism, then 
we apply it to the WMAP 5-year temperature map de- 
scribing the features we observed and their significance 
in Sec. III. In Sec. IV we address the angular power 
spectrum modification induced by the specific features 
we measure in the temperature map, computing the ef- 
fect on the cosmological parameters. We then draw our 
conclusion in Sec. V. 



II. NEEDLETS AND SEARCH FOR FEATURES 

We provide here a brief summary of the needlet con- 
struction, while we refer to [26] and the references therein 
for a more complete discussion. We recall the spherical 
needlet function is defined as 



where {^j/e} are the cubature points on the sphere, corre- 
sponding to the frequency j and the location /c, and 6(.) is 
a filter function in harmonic space. Each cubature point 
requires a proper weight, Aj^. Since our implementation 
is based on the Healpix pixelization^ [40], we actually 
compute Eq. 1 on the center of the pixels for a given 
resolution, following the procedure described in [25]. In 
this context the cubature weights can be approximated 
by 1/A/p, where is the number of pixels for the chosen 
Healpix resolution and k corresponds to a Healpix pixel 
number. 

In general, needlets can be viewed as a convolution of 
the projection operator Y.rn^ irninY^irn^ijk) with suit- 
ably chosen weights provided by the function 6(.). Details 
on the function 6(.) we used can be found in Ref. [26]. 
We show its profile in Fig. 1 for several values of in 

the sequel we will use the notation hi iox hi^-^^. B is 

a user-chosen parameter that characterizes the weights 
b£ and then the entire set of needlets, since its value de- 
termines the width of the filter function. The choice of 
B must be driven by the insight on the range of mul- 
tipoles to be probed. Since we are interested in large 
angular scales, we set B = 1.8: this choice allows us to 
have eight frequencies spanning multipoles up to £ = 200, 
while the information at smallest scales is concentrated 
on just three frequencies. The range of multipoles cov- 
ered by each needlet is summarized in Table I. 

Since the function has a finite support in £-space, 
by construction each needlet receives a contribution only 
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FIG. 1: Top: Profile of the function b(x) in ^-space for the 
choice B = 1.8. The black solid line represents j = 4, the 
red dot-dashed Une j = 3 and the blue dot-long-dashed Une 
j = 5. Bottom: needlets on the sphere for j — 4. At each 
cubature point, (fjfc, the needlet is sharply localized. 



from a specific range of multipoles, which is related to 
the index j by a logarithmic function oi B. In pixel 
space needlets are quasiexponentially localized, due to 
peculiar properties of the window function 6^, as shown 
by [23] . We are thus able to exploit a tight control on the 
localization both in terms of angular scales and directions 
on the sky, thus keeping a record of the signal in both 
domains. This makes possible the detection of specific 
features in the CMB map that are contributing power, 
for example, over a well-defined range of multipoles of 
the angular power spectrum. 

To this aim, we first decompose the CMB temperature 
map onto the needlet system. As recalled before, needlets 
do not make up a basis, but a tight frame; however the 
latter enjoys all the reconstruction properties that are 
usually associated to bases (indeed tight frames can be 
simply viewed as a redundant basis). The needlet coeffi- 
cients Pjk result from the projection of the temperature 



field on the elements of the needlet system: 

= ^^rnYimi^jk)^ (2) 



j 


range 


1 


2 


2 


2-5 


3 


4-10 


4 


6-18 


5 


11-33 


6 


20-60 


7 


35-108 


8 


63-196 


9 


113-352 


10 


203-635 


11 


365-1143 



TABLE I: Range of multipoles spanned by needlets for B — 
1.8. 



III. MAPS 

We decomposed the Internal Linear Combination 
(ILC) WMAP 5-year temperature map by means of 
needlet functions corresponding to 5 = 1.8. As discussed 
in the previous section, this specific choice provides us 
needlets at 11 frequencies j, that span properly the low 
multipoles we are interested in. As described above, the 
expression 2 can be easily represented by a map in the 
Healpix scheme, in which each pixel is associated with 
the needlet coefficient evaluated at ^j/^. Figure 2 shows 
three clarifying examples of the WMAP 5-year ILC map 
decomposition on j = 3, j = 4 and j = 5 needlets by 
making use of the functions bi depicted in Fig. 1. 

We first applied our procedure to the WMAP map with 
the extended mask KQ75 from the WMAP team. Such 
a large mask, however, entails a greater correlation for 
multipoles corresponding to large angular scales; for low 
values of j, this could potentially impact the needlets' co- 
efficients as well. Any detection of significant features in 
the CMB anisotropy map using the large mask, however, 
can be reconfirmed at a higher significance level with a 
smaller mask. After studying the temperature map using 
KQ75, we repeated all the analysis applying KQ85, which 
is the mask currently favored for cosmological data anal- 
ysis, including power spectrum measurements. We find 
that our results related to statistically significant features 
on the WMAP map are fully consistent in the two cases. 
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III. A Hot/Cold spot maps 

In Fig. 2, we represent the needlet coefficients for the 
cases of j = 3 to j = 5. These three frequencies are par- 
ticularly interesting since they probe the low-multipole 
region, highlighting a peculiar pattern of anisotropics in 
the southern hemisphere. When j = 4, we recognize 
two spots, one hotter and one colder than the average 
CMB fluctuations. While the latter is a detection of the 
so-called cold spot^ a feature studied in detail in the lit- 
erature as a source of non-Gaussianity in the CMB map 
[5, 41-46], the hot spot has not had the same scrutiny 
with only a minor description in Refs. [43, 47]. 

In fact, the j = 3 and j = 5 needlet coefficients 
represent a cross-validation of the existing literature, 
as they also highlight a second hot spot centered at 
{gi = 173,^5 = —46) (as pointed out in Ref. [47]) and 
a minor cold spot centered at {gi = 80,^5 = —33) (ob- 
served in Ref. [43]), that actually appears with an adja- 
cent hot spot. In Fig. 3 the three spots we used in the 
angular power spectrum analysis are shown. We stress 
that our identification follows uniform criteria with the 
same technique, quite differently from some of the exist- 
ing literature. Moreover, as explained above we are also 
able to identify exactly the range of multipoles where 
these features contribute to the angular power spectrum 
of CMB anisotropics, due to the specific needlet proper- 
ties. 

To establish the significance of the features we de- 
tected, we consider a Monte Carlo analysis by performing 
a large set of simulations (1000) of a Gaussian CMB sky 
with the same angular power spectrum as the WMAP 5- 
year best-fit cosmological parameters [48] ; we then com- 
pute the average and the standard deviation of the dis- 
tribution for each needlet coefficient. The expected dis- 
tribution is Gaussian with zero mean, (see [24]), in good 
agreement with our simulations. We focus on the statis- 
tic 



\Pjk- < Pjk > I 



(3) 



where ap.,^ is the usual standard deviation of the distri- 
bution. In Fig. 2, we show that the hot and cold spots 
exceed three sigma level at j = 4. The cold spot appears 
significantly both in the needlet coefficients for j = 4 and 
j = 5, and thus its impact can be reckoned to span the 
range between i = 6 and i = 33 (Table 1). In Table II 
the significance values of the three anomalous spots are 
quoted. 



Property / Spots 


Cold Spot 


Hot Spot 


Hot Spot 




(209, -57) 


(320, -34) 


(173, -46) 


j 


4 


4 


3 


nmax 


(-)3.72 


3.56 


3.24 



TABLE II: Main properties of the spots highlighted in our 
analysis. 



It may be suspected that the hot spot we located could 
be a spurious effect due to oscillations in the needlet func- 
tion. We considered this explanation, but we concluded 
that the distance at which the hot and cold spot appear 
is greater than the needlet oscillation range. We aim at 
a further investigation of this issue in our ongoing work. 

As a test for the joint significance of needlet coeffi- 
cients, we consider the statistic 



^ k 



PjkPj'k- < PjkPj'k > 



^k 



(4) 



where E-^-^ is the second moment of the distribution. Of 
course, the statistical distribution of F cannot be taken 
as Gaussian. Hence, to analyze its statistical significance 
we used again a set of Monte Carlo simulations and we 
followed the procedure described in Sec. III. Figure 4 
shows the values of these statistics for the pair j = 3 and 
j = 4, left panel, and for j = 4 and j = 5, right panel. 
Results are expressed directly in terms of the statistic 
defined in Eq. 4 in the upper panels; in the lower panels 
a threshold of 7.5 for F;^-^ is adopted to underline the fact 
that the hot and cold spots identified by the mask j = 4 
are the most significant. 



III.B North-South Asymmetry 



We now wish to investigate the extent in which mask- 
ing the previously found features affects the asymmetry 
between the northern and the southern hemisphere of the 
CMB sky. To see this, we recall that in Refs. [24-26] it 
is shown that 



47r 



whence jSj = is an unbiased estimator for the 

weighted angular power spectrum. In [24], Sec. 7, it is 
also shown that this statistic is approximately Gaussian 
(after centering and normalization) at high frequencies, 
see also [30, 32] for extensions and related work. Devel- 
oping this idea, we computed here the total power in each 
hemisphere by taking the sum of squares of the needlet 
coefficients at j = 3 and j = 4 extracted from the masked 
ILC temperature map 



(5) 



J=3 k 



J=3 



We label the mask associated with the j = 3 and j = 4 
features we found as "j3j4" , and we include it in addition 
to standard WMAP masks. 
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FIG. 2: Needlet coefficients of tfie WMAP 5-year CMB temperature map. On tfie left, from top j = 3, j = 4 and j — 5 are 
plotted. The set of needlets is characterized hy B = 1.8. Each pixel represents the coefficient for the needlet function computed 
at ^jk, where k identifies the pixel in the Healpix ordering. The effect of the applied KQ75 mask is clearly visible, setting to 
zero the value of each pixel that belongs to the mask. It is interesting to notice that needlet coefficients highlight the presence 
of the well-known cold spot in the southern hemisphere, as well as a hot spot localized in the southern hemisphere closer to the 
mask. Needlet coefficients for j = 5 show the cold spot pretty well, while the hot spot is weaker. Another couple of hot/cold 
spots appear above the Galactic plane. On the right, from the top, significance of the needlet coefficients for j = 3, j = 4 and 
j = 5. The three maps show the significance Sjk above the threshold of 3. This allows us to localize in ^-space the contribution 
of the hot spot that results to be in the range of multipoles between = 6 and ^ = 18. Computing the coefficients for j = 6 
and observing that the cold spot, if present, does not have a high statistical significance, we can deduce the range of multipoles 
covered by the cold spot: between i = 6 and ^ = 33. 



The difference is measured by computing 

D^^ (C|4 - C^^a) . (6) 

where C3 4, i = iS*, TV is the quantity in Eq. 5, where the 
needlet coefficients are restricted to either the northern 
or the southern hemisphere (in the Galactic coordinate 



system), normalized to the variance V of the whole sky. 
The latter actually turn out to be the cosmic variance 
of the CMB power spectrum binned with bj In some 
sense, D is measuring the difference between two local 
versions of the power spectrum estimator; such local es- 
timators can indeed be rigorously justified, in view of the 
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Hot Spot 
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FIG. 3: Detail of the spots. From the top, the cold spot at 
{gi = 209, Qb = -57), the hot spot at (gi = 320,^6 = -34) 
present at j = 4, and the second hot one at {gi = 173, gt = 
—46) measure at j = 3. The true colors are altered by the 
use of the masks we employed in the analysis to highlight the 
region above three sigma level. 



to the northern one that is reduced by a factor of 2 when 
the "j3j4" mask is applied. Notice that, as expected, the 
power measured in the north region is not affected by 
the masking procedure. In the lower panel of Fig. 5 we 
quantify this effect computing the difference — 
normalized to the variance of the estimator. 

The results are summarized in Table III. Our findings 
are sensitive to the chosen sky-cut: when the KQ85 mask 
is used, the asymmetry, measured in terms of a difference 
in the variance of power, decreases from 4.26 without the 
j3j4 needlet mask to 1.97 with needlet mask. When the 
more aggressive mask KQ75 is applied to CMB data (to- 
gether with the j3j4 mask), the difference is larger and 
the global asymmetry is further reduced: from 4.37 to 
2.0. Note that the north-south power variance difference 
with just the WMAP team's masks, KQ75 and KQ85, is 
rather small (4.26 vs 4.37); this suggests that the hemi- 
spherical asymmetry cannot be explained by simply ex- 
tending the galactic plane mask. With the j3j4 mask we 
introduced (which cuts roughly 0.5% of the sky), we find 
a significant reduction of a factor of 2. 

Given that we mask a smaller area on the sky than 
the 10% difference between KQ85 and KQ75, it seems 
rather likely that the hemispherical asymmetry can be 
credited to features in the southern hemisphere. While 
we have localized (some of) these features, this does not 
establish by itself whether the asymmetry is primordial or 
associated with fluctuations in our local universe. There 
exist already several extended studies on the nature of 
the cold spot ([45, 49-54]); for further statistical studies 
we make publicly available the j3j4 mask. 



Mask 


^3,4 


^3,4 


^3,4 


D 


kq85 


766 


556 


974 


4.26 


kq85-hj3j4 


673 


544 


802 


1.97 


kq75 


703 


502 


908 


4.37 


kq75-hj3j4 


655 


492 


728 


2.0 



TABLE III: We report the values of total power (Eq. 5) car- 
ried by needlets at j = 3 and j — ^ extracted on the ILG 
map. Four cases, corresponding to the different masks we ap- 
plied in this analysis, are shown. The last column reports the 
significance D as defined in Eq. 6. It is interesting to notice 
how masking the hot and cold spots reduces the asymmetry 
by a factor greater than 2, while the 'y3j4" mask covers 0.5% 
of the sky only. 



uncorrelation properties of needlets in pixel space (see 
[28]). In practice, we defined the pixels in the northern 
hemisphere as those outside the mask characterized by 
6 < tt/2 and pixels in the southern hemisphere as those 
outside the mask, with 9 > 7r/2. 

In Fig. 5 we plot (3j as defined in Eq. 5 extracted from 
the whole sky as well as those measured in each hemi- 
sphere, both including the spots and masking them. The 
southern hemisphere shows an excess of power compared 



III.C Further statistical analysis 

As a further statistical cross-check, we chose randomly 
needlet coefficients at different locations on the map and 
fit their distribution to the one derived from 1000 sim- 
ulations where the KQ85 mask was applied. The simu- 
lated results are in excellent agreement with a zero-mean 
Gaussian distribution, as shown in Fig. 6. This result is 
of course expected, as the needlet coefficients are a linear 



FIG. 4: r statistics for j = 3 and j' = 4 (left-hand side) and with j = 4 and / = 5, (right-hand side). Both the cold spot and 
the hot spot are clearly visible in the maps, but when a threshold of Tj,^ > 7.5 is applied their signal is strong only on the 
j = 4 and / = 5 maps. 



functional of the underlying temperature map. However, 
we report the figure as a further check to assure that the 
procedure we followed to compute the significance of the 
spots is well justified. 



We followed the same procedure also to quantify the 
significance of Tj^^ . Of course, in this case simulations 

are indeed necessary, because P-^-^ is a nonlinear statistic 
and hence non-Gaussian. In Fig. 7 we provide some ev- 
idence on the significance of the statistics we measured 
in the regions where the anomalous spots are located. 
The curve is the fit to the distribution of needlet coef- 
ficients in simulated maps, while the vertical lines mark 
the value measured in the WMAP temperature needlet 
coefficients. The simulated distribution has rather large 
tails; the statistical significance of F-^-^ at the locations 
corresponding to the spots is nevertheless rather high, 
with an estimated p- value of 0.5% for the cold spot and 
of 1.65% for the hot spot. We can still confirm that the 
signal is mostly peaked in the correlation between j = 4 
and j = 5. 



IV. IMPACT ON THE CMB POWER 
SPECTRUM 

The purpose of this section is to investigate the ex- 
tent in which masking the hot and cold regions we found 
affects the GMB power spectrum [49]. To address this 
issue, we first estimate the angular power spectrum from 
the ILC WMAP 5-year map after applying the KQ85 
mask. We then compare this result with the angular 
power spectrum resulting from a wider mask: the sum 
of KQ85 plus the regions above three sigma level we dis- 
covered when performing our temperature analysis. We 
report in Fig. 8 the resulting mask. 

The effect of the different masking is not negligible, 
reaching the value of 12% at low multipoles. To check 
against systematics, we performed a Monte Carlo sim- 
ulation of 200 CMB maps with the underlying WMAP 
5-year best-fit model [48] , in order to estimate the mean 
effect of the joint mask. We computed the average and 
the standard deviation to quantify the hot/cold spots ef- 
fect. The results are shown in Fig. 9. 

Quite remarkably, the region where the signal is 
stronger is exactly the one where the cold spot is local- 
ized. This may be interpreted as a confirmation of the 
localization properties of needlets in pixel and harmonic 
spaces. We believe the use of the ILC map is justified 
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FIG. 5: In the upper panel pj extracted from WMAP data 
when different sky cuts are apphed are shown. The black solid 
Une shows the signal for the whole CMB sky (kq85 is applied) . 
The blue and red dashed lines show how the power is split 
between the two hemispheres. When the cold and hot spots 
detected are masked the excess of power in the southern region 
is decreased (light blue and orange dot-dashed lines). In the 
lower panel the significance of the difference D is plotted. 



here, because our signal peaks at low multipoles; how- 
ever for completeness we computed the same quantity 
from both the W and V bands of WMAP 5 year, as well 
as using the map obtained by [55]^. The result is fully 
consistent with what we found using ILC, thus validating 
the procedure we followed. The signal is shown in Fig. 10 

The next step has been to evaluate the effect that the 
change in the angular power spectrum induces on cos- 
mological parameter estimates. Since the low multipoles 
region is affected by the largest variation, we expect that 
changes might occur on the spectral index, n^, the optical 
depth, r and possibly on the primordial fluctuation nor- 
malization amplitude Ag. Actually, since the variation is 
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FIG. 6: From top distribution of the needlet coefficient for 
J = 3, J = 4 and j — 5. The blue (on the left-hand side) and 
red (on the right-hand side) vertical lines mark the values 
of the cold and hot spot respectively. They are both well 
beyond three sigma level for j = 4. At j = 3 he significance is 
basically zero, while at j = 5 only the cold spot is still visible. 



at maximum 12% in an handful of multipoles, we expect 
a global variation on the power spectrum of roughly few 
parts on 1000. 

The WMAP team performed the temperature analysis 
splitting the low multipoles and the high multipoles re- 
gions. The former is probed by a Gibbs-sampling based 
Monte Carlo analysis [56] . The high moments are instead 
investigated fitting the angular power spectrum extracted 
from the W and V bands. To take into account the mod- 
ifications due to the new masking we replaced the KQ85 
mask used at low resolution with the joint mask KQ85 
plus the hot and cold spots (Fig. 8), and multiplied the 
angular power spectrum used for the analysis by the ra- 
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FIG. 7: Distribution of T statistic for (j,/) = (3,4) (top) and 
— (4, 5). The two vertical lines mark the values of the 
hot and cold spot, the latter being the more significant. As 
expected, this is non- Gaussian, and is characterized by large 
non-Gaussian tails. 



KQ85 + Hot/Cold Spot mask 
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FIG. 8: Joint KQ85 and hot /cold spots mask applied to the 
WMAP ILC temperature map. 



tio between the Ci computed with the wider mask and 
Ci obtained with the unmodified KQ85 mask applied to 
ILC, in the range of multipoles 2 — 200 (see Fig. 10). 

The results are shown in Table IV. In short, we do not 
observe a significant variation in any of the cosmological 
parameters. 

It is well known that the angular power spectrum mea- 
sured by the WMAP Collaboration shows some interest- 
ing features at low multipoles; in particular the range 



Effect of HOT/COLD spots 
1.2 F ' ' ' ^ 




FIG. 9: Effect on the angular power spectrum due to the 
subtraction of the hot and cold spots in the CMB temperature 
map. The red solid line shows the difference in the Ci] the 
gray dashed-dotted line represents the average modification of 
the simulation, while the solid lines mark the one sigma level. 
In the region between i = S and i = 30 the effect exceeds one 
sigma level: that region is the one where the effect of the cold 
spot is stronger. 



Mask effect: map comparison ~ kq85 




FIG. 10: Variation due to the combination of the hot and cold 
spots on five different power spectra: ILC the yellow curve, 
W the black curve, V the blue one, WV the light blue curve, 
and that extracted by the map reconstructed in [55] in orange. 
Except for the WMAP ILC map at ^ > 200, the difference 
on the angular power spectrum is consistent among a variety 
of maps. Beyond i = 200 the ILC power spectrum shows 
features not compatible with other spectra, probably due to 
the way the different WMAP channels are combined. 



between ^ = 20 and £ = 24 has a deficit in power with re- 
spect to the prediction of the best-fit ACDM theoretical 
model, while that between i = 37 to i = 44 shows excess 
power. To investigate these issues, we chose needlets cor- 
responding to frequencies that match those two intervals 
and we looked for coefficients exceeding the threshold of 
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Parameter 


WMAP5 


Hot/Coid spot masked (j4j 


J 3- J 4 mask 




0.0227 ± 0.0006 


0.0228 ± 0.0006 


0.0228 ± 0.0006 




0.110 =b 0.006 


0.109 ± 0.006 


0.109 ± 0.006 


Oa 


1.040 ±0.003 


1.040 ±0.003 


1.040 ±0.003 


T 


0.089 ±0.018 


0.091 ±0.017 


0.089 ±0.017 


ns 


0.965 ±0.014 


0.966 ±0.014 


0.966 ±0.014 


In(lO^^As) 


3.18 ±0.05 


3.17 ±0.05 


3.17 ±0.05 



TABLE IV: Effect of the wider mask on the ACDM six parameters. The difference due to the sum of KQ85 mask plus hot/cold 
spot mask is small. 



three sigma. More precisely, we select B = 1.2 and we 
take j = 17, j = 20 to span the relevant ranges of multi- 
poles. Figure 11 (upper panel) and Fig. 12 show (respec- 
tively) the b£ profile employed for this purpose and the 
corresponding needlet coefficients. 

Profile ~ B = 1.2 
1.2 r ' ' ' ' 1 



0.8 - I 1 : : 

0.6 - I i ! ': 

0.4 - I ■ ! 

0.2 - ■ i I ! 

0.0 r I ' , \ 

1 10 100 

1 




kq85 W band 
kq85+j W band 



1 10 100 



FIG. 11: Top panel, profile of the function b{x) in ^-space 
for the choice B — 1.2. The red dot-dashed line represents 
j = 17 and the blue long dashed line j = 20. Lower panel, 
power spectrum modification due to the structures measured 
using the set of needlets shown in Fig. 12. 

We find that the effect on the angular power spectrum 
is actually rather small. However it does follow the ex- 
pected sign, decreasing the deficit and the excess of power 



in the selected ranges. The effect is shown in the lower 
panel of Fig. 11. 



V. CONCLUSIONS 

We apply spherical needlets to the Wilkinson Mi- 
crowave Anisotropy Probe 5-year cosmic microwave back- 
ground dataset, to search for imprints of nonisotropic 
features in the CMB sky. 

After calibration by means of a large set of mock simu- 
lations, the analysis of needlet coefficients highlights the 
presence of the now well-known "cold spot" of the CMB 
map in the southern hemisphere, and in addition two hot 
spots at significance greater than 99% confidence level, 
again in the southern hemisphere and closer to the Galac- 
tic plane. While the cold spot primarily contributes to 
the anisotropy power spectrum in the multipoles between 
^ = 6 to = 33, the hot spots are found to be dominat- 
ing the anisotropy power in the range between £ = 6 and 
i= 18. 

We also studied the effect the two spots have on the 
CMB power spectrum, by building 1000 mock CMB sim- 
ulations. We conclude that, especially at low multipoles, 
the effect is measurable: masking both the cold and the 
two hot spots results in an increase in the quadrupole am- 
plitude of 10%, while at ^ = 10 power is reduced by 12%. 
To investigate the effect of this difference on the value of 
cosmological parameters, we modified the WMAP 5-year 
CMB fiducial power spectrum and the KQ85 mask used 
to perform the analysis by cutting out the contribution 
of the two spots, and we repeated the parameter estima- 
tion analysis. The results we obtain are slightly different, 
but fully consistent within the la errors on parameters 
published by the WMAP team. 

Since all three spots appear in the southern hemi- 
sphere, we also studied the power spectrum asymmetry 
between the two hemispheres, which has been previously 
found to be statistically significant. When the features 
detected by needlets are masked, we find that the differ- 
ence in the power, measured in terms of the anisotropy 
variance between £ = 4 and £ = 18, is reduced by a 
factor of 2. This decreases the significance of the pre- 
viously claimed north-south asymmetry. We make the 
mask resulting from needlet features available for future, 
more detailed studies on the asymmetries in the CMB 
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WMAP 5yr Temperature map Necdlets coefficients Significance: j = 17 
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FIG. 12: Needlet coefficients and their significance at j = 17 and j = 20 for B = 1.2. 
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